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Kink ratchets in the Klein-Gordon lattice free of the Peierls-Nabarro potential 
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A discrete Klein-Gordon model with asymmetric potential that supports kinks free of the Peierls- 
Nabarro potential (PNp) is constructed. Undamped ratchet of kinks under harmonic AC driving 
force is investigated in this model numerically and contrasted with the kink ratchets in the conven- 
tional discrete model where kinks experience the PNp. We show that the PNp-free kinks exhibit 
ratchet dynamics very much different from that reported for the conventional lattice kinks which 
experience PNp [see, e.g., Phys. Rev. E 73, 066621 (2006)]. Particularly, we could not observe any 
significant influence of the discreteness parameter on the acceleration of PNp-free kinks induced by 
the AC driving. 

PACS numbers: 05.45.Yv, 63.20.Ry, 05.60.Cd 



I. INTRODUCTION 

Ratchet dynamics of a point-like particle or a quasi- 
particle such as soliton is the motion of the particle in a 
certain direction under AC force whose average is zero. 
Ratchet transport phenomenon can be observed under 
the following two conditions: (i) the system must be out 
of thermal equilibrium and (ii) the space-time symme- 
try of the system must be violated [IL B 0] ■ Ratchet 
dynamics has been receiving growing attention of re- 
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state physics |18| . 

Soliton ratchets were first studied by Marchesoni [l^ 
for the overdamped Klein-Gordon equation. Unlike 
point-like particles, solitons can have internal vibra- 
tional modes [l^] and these modes can strongly affect 
the ratchet dynamics especially for underdamped case 
pH . [2^ . So far, ratchet dynamics have been stud- 
ied for a number of continuum soliton bearing sys- 
tems [n nn, m m [13, ii m, Hi Si, whue in 

many applications ratchet was observed for discrete kinks 
[ia, 14, ij, p. Effects of discreteness on kink ratch- 
ets have been studied by Zolotaryuk and Salerno [29j . 
They have found that in comparison with the continuum 
case, the discrete case shows a number of new features: 
nonzero depinning threshold for the driving amplitude, 
locking to the rational fractions of the driving frequency, 
and diffusive ratchet motion in the case of weak inter- 
site coupling. For the damped, driven Frenkel-Kontorova 
chain, which is the discrete analog of the sine-Gordon 
equation, Martinez and Chacon have shown that phase 
disorder introduced into the asymmetric periodic driving 
force can enhance the ratchet effect [30]. 

In discrete systems translational invariance is typically 
lost and static solitary waves usually cannot be placed 
arbitrarily with respect to the lattice but only in the 
positions corresponding to the extremums of the Peirls- 
Nabarro potential (PNp), induced by the lattice. Con- 



figurations corresponding to the maximums of PNp are 
unstable while those corresponding to the minimums are 
stable. It has been found that the presence of PNp makes 
kink ratchets much more complicated in comparison to 
the continuum case [29]. 

On the other hand, in the recent past, several differ- 
ent families of discrete Klein-Gordon systems supporting 
translationally invariant static solutions with arbitrary 
shift along the lattice have been derived and inve stig ated 

P, ii, IM H ii, [11, il [11, [M EO, IM El, 

Such discrete models are often called exceptional. Phys- 
ical properties of solitary waves in the exceptional dis- 
crete models were found to be very much different from 
their conventional counterparts. For instance, they can 
support conservation of momentum 3l|, 38 and can sup- 
port kinks which move with a special }44| or arbitrary 
|45| velocity without emitting radiation. Peculiarities of 
kink collisions in such models have been investigated in 
the work (ssj . Static solutions in the exceptional dis- 
crete systems possess the translational Goldstone mode 
(3^ [35I . This means that such static solutions are not 
trapped by the lattice and can be accelerated by arbi- 
trary weak external field. Exceptional discrete models 
can describe physically meaningful systems [4l[, that is 
why investigation of physical properties of such systems 
is very important. 

In the present study we continue the investigation of 
physical properties of the exceptional discrete models 
that support static kinks free of PNp. The main goal 
of the study is to see how the special properties of the 
kinks can influence their undamped ratchet dynamics un- 
der single-harmonic driving. 

In order to achieve this goal we need to construct 
an exceptional discrete Klein-Gordon model with asym- 
metric background potential. It is important that the 
constructed model be Hamiltonian, otherwise energy in- 
crease in the system can be observed even for small- 
amplitude driving with the frequency laying outside the 
phonon band. 

The paper is organized in five Sections. In Sec. [TTl 
we describe two discrete, Hamiltonian Klein-Gordon sys- 
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terns, with and without PNp, having asymmetric on-site 
potential and having common continuum hmit. In Sec. 
mil we compare the properties of static kinks in these two 
models and then in Sec. IIVI study kink ratchets mostly 
for the PNp-free model. Section IVl concludes the paper. 

II. DISCRETE KLEIN-GORDON MODELS 
WITH ASYMMETRIC POTENTIAL 

A. General formulation 

The Klein-Gordon field has the Hamiltonian 




where 4>{x, t) is the unknown field and V{(j)) is a given po- 
tential function. The corresponding equation of motion 
is 

(t>tt = <j^.x - V'{c^) , (2) 

where V'{(j)) = dV/dc/). 

Equation ^ will be discretized on the lattice x = nh, 
where n — 0, ±1. ±2..., and h is the lattice spacing. Tra- 
ditional discretization of Eq. ([2]) is 

<i>n = ■J^{4'n-1 - 2(/)„ -I- 0,1+1 ) - V'{(l)n)- (3) 

Using the discretized first integral (DFI) approach offered 
in [SJI and developed in one can construct a discrete 
model whose static version is an integrable map. Follow- 
ing this method we begin with the first integral of static 
Eq. (El), 0^ - 2F(0) + C = 0, where C is the integra- 
tion constant. The first integral can also be taken in the 
following modified form [33| 

v{x) = 0, ~ y/2VW~C = . (4) 

Next step is to rewrite the Hamiltonian Eq. ^ in terms 
of v{x) as follows 

H=2j l"^* + + 2<jy,^2V{cj>)~c]dx , (5) 

where we omitted the constant term. 

The first integral Eq. ([4]) can be discretized as follows 



v{4>n-lAn) = ^" ' - ^2V{^n~u4^n)-C = 0,(6) 

where we demand that V{4>n~i, <t>n) — * y{4>) in the con- 
tinuum limit {h — > 0). Thus we obtain the discrete ver- 
sion of the Hamiltonian Eq. ([T]) 

11 

+ 2 '^"~^'^"-V 2n'^n-l>n)-C}. (7) 



Final step is to discretize the background potential as 
suggested by Speight [sll , 

/— — — G(0n) - G(0n-l) 

V2l/ (0„-l,0„) - C = , 

(Pn - (pn-l 

where G' (0) = ^2^(0) - C. (8) 

With this choice the last term of the Hamiltonian Eq. ([7]) 
reduces to (2//i)[G((?!)„) — G'((/)„_i)] and it disappears in 
the telescopic summation. Further, according to Eq. ([8]), 
the discretized first integral Eq. ^ assumes the form 

/ , , N 0n - G((/>„) - G{(t)n-l) , . 

W(0„_1,0„) = , (9) 

" (Pn- 0n-l 

and the equations of motion derived from Eq. ([7|) with 

v{cl)n-lAn) ^ v{(f)n^i,(j>n) are 

d 

(pn = -v{4>n-l,4'n)Tr—v{4>n-l,4'n) 
0(pn 

d 

-v{(j)n,(j)n+l)T-—v{(j)n,4>n+l)- (10) 

Obviously, equilibrium static solutions of this model 
can be found from the two-point nonlinear map 
v{(j)n-i,4>n) = 0, where v{4>n-i, 4>n) is given by Eq. ([9]). 
Such solutions can be constructed iteratively starting 
from any admissible initial value (pn-i or and thus, 
the PNp is absent for such family of equilibrium solu- 
tions. 

B. Polynomial asymmetric potential 

Klein-Gordon kink ratchets are possible if the on-site 
potential or the driving force or both are asymmetric. We 
study the kink ratchets under single-harmonic AC driv- 
ing and thus, the on-site potential must be asymmetric. 

We take G' in Eq. ([8]) in the form of the quartic poly- 
nomial function 

G'{(f>) ^a(f>'^ + b(f)^ + c(f> + d, (II) 

where the cubic term was not taken into account because 
it can always be removed by a proper shift cj) cj) ~ 
Then the on-site potential (with C = 0) is 

V{^) = i (a(/.4 + ^ c0 + d) ' . (12) 

The simplest discrete Klein-Gordon model corresponding 
to this potential (will be referred to as DKGMI), accord- 
ing to Eq. ©, is 

4>n = (0n-l - 2(/)„ + (pn+l) 

- (a0^ + hcf>l + c0„ + d) {Aa4>l + 260n + c) , (13) 
and its Hamiltonian is 

n 

+ (a(/.^ + 602 +C(/.„ + d)'}. (14) 
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FIG. 1: On-site potential Eq. for 6 = and a = ~1, 

c = 0, d = 1 (dashed line); a = -0.83988, c = 0.382925, d = 
0.860756 (dash-dotted line); a = -0.643049, c = 0.626843, 
d = 0.706344 (solid line). 

A more sophisticated discrete Klein- Gordon model 
(will be referred to as DKGM2) is defined by Eq. ITOl) 
with 




b c 

-3 (<^n + Mn~l + 4>l-l) - 2 (-^^ + -^"-l) - '^^ (15) 

which was found by substituting Eq. (|lip after integrat- 
ing into Eq. This model has the Hamiltonian 

^2 = ^E{'^" + K0«-i,0„)f }. (16) 

n 

The asymmetry of the potential Eq. is controlled 
by the parameter c and for c = the potential is symmet- 
ric. The on-site potential Eq. has four parameters. 
Let us fix the height of the potential barrier equal to 0.5, 
the distance between the two minima equal to 2. The 
asymmetry can be chosen by setting a value of c. There 
is still one free parameter and we will take & = 0. 

In Fig. [T] we plot the on-site potential Eq. \12\ for 
6 = and a — —1, c = d — \ (dashed line); a = 
-0.83988, c = 0.382925, d = 0.860756 (dash-dotted line); 
a = -0.643049, c = 0.626843, d = 0.706344 (solid line). 
One can see that the asymmetry of the potential increases 
with c. 

The asymmetric potential supports two vacuum so- 
lutions, 0„ = 4>i and (pn — 4>2, where 0i and 02 are 
the coordinates of the two minima of the on-site po- 
tential. Small-amplitude phonon vibrations of the form 
4>n ~ exp[i(qn — ujt)], where q is the phonon wavenumber 
and uj is the phonon frequency, have different spectra for 
different vacuums. 

Borders of the phonon bands for each vacuum can be 
found for DKGM2 from 



ujI = 28a20f + 30ab(j)j + 2Oac0^ 
-^3 (26^ + Aad) (Pi + &hc<i), + + 2M, (17) 

28 4 

2 \ 2 2 4 

+6 f y + gfldj 4\ + 3M. + 3^^^ + (18) 

where j — 1^2 and (aj2) corresponds to g = (g = tt). 

For DKGMl the borders corresponding to q = co- 
incide with that for DKGM2, while the borders corre- 
sponding to g = TT are 

ul = 28a20^ -I- 3Oa60^ -I- 2Oac0^ + 

6 (6^ + 2ad) (Pj + &bccj), + + 2bd+ ^. (19) 

Numerical results in this work will be obtained for the 
on-site potential with the parameters a = —0.643049, b = 
0, c = 0.626843, d = 0.706344 (shown by the solid hne 
in Fig. [1]). The potential has minima at = —0.768678 
and 02 — 1.231321 and a maximum at 0max = 0.62462. 



III. PROPERTIES OF STATIC KINKS IN TWO 
LATTICES 

Before we proceed with the we need to study the prop- 
erties of the kinks in DKGMl and DKGM2 because they 
will help us to interpret the results of the kink ratchet 
dynamics study. 

Equilibrium static kink solutions for the DKGMl can 
be found numerically while for DKGM2 they can be 
found iteratively using Eq. (jl5p for any initial value of 
4>n (or (j)n-i) lying between two minima of the on-site 
potential. 

In the DKGMl there exists only one stable static kink 
configuration [shown in Fig. [2] (a)], corresponding to the 
minimum of the PNp. Static kinks in the DKGM2 do 
not experience the PNp and they can be placed anywhere 
with respect to the lattice. A family of equilibrium kinks 
is presented in Fig. [2] (b). Kinks in both models have 
asymmetric shape because of the asymmetry of the on- 
site potential. 

Small-amplitude oscillation spectra for the chains con- 
taining a kink are presented in Fig. [3] for different values 
of the discreteness parameter h. Dashed horizontal line 
shows the lower edge of the phonon band which is the 
same for both DKGMl and DKGM2 and can be found 
from Eq. (|17p for J = 1 (soft minimum) . Presented spec- 
tra show the kink's internal vibrational modes. Since 
both discrete models share the same continuum limit [de- 
fined by Eq. ([2]) and Eq. p2)) ]. for small discreteness 
{h < 0.25) their spectra are close. Kink in DKGMl [see 
in (a)] possesses two internal modes, one of them is the 
destroyed translational mode (for small h it approaches 
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FIG. 2: (a) Stable equilibrium kink in DKGMl and (b) a 
family of equilibrium kinks in DKGM2. In both cases h = 
0.6. Dashed lines show the locations of minima of the on- 
site potential and the dash-dotted line the location of the 
maximum. 





h h 



FIG. 3: Kink's internal mode frequencies as the func- 
tions of the discreteness parameter for (a) DKGMl and (b) 
DKGM2. Dashed line shows the bottom edge of the phonon 
band. Model parameters (here and in the following); a — 
-0.643049, 6 = 0, c = 0.626843, d = 0.706344. 



zero frequency). Kink in DKGM2 [see in (b)] possesses 
the zero-frequency translational mode for any /i, and for 
h > 0.48 two new internal modes appear. Note that 
the spectrum in (b) was calculated for the kink having a 
particle at the maximum of the on-site potential and the 
kink internal mode frequencies within the studied range 
of h are only shghtly dependent on the location of the 
kink with respect to the lattice. For example, maximal 
difference between the internal mode frequencies for dif- 
ferent kink positions with respect to the lattice observed 
at /i = 1 is 0.9%. 



FIG. 4: (a) Kink motion in DKGM2 at h = 0.6 for two dif- 
ferent values of initial phase of driving force if at the same 
amplitude A — 0.04 and frequency Q — 0.5. Oscillating lines, 
presenting the kink coordinate as the functions of time, are 
fitted by the square parabolas, (b) Acceleration of the kink 
and (c) initial velocity of the kink as the functions of ip. 

IV. KINK RATCHETS 

To study kink ratchets we add to the right-hand sides 
of equations of motion Eq. ^ (DKGMl) and Eq. ^ 
(DKGM2) the harmonic external force 

F{t) = Acos{nt + ip), (20) 

with the amplitude A, frequency Q, and initial phase ip. 

The initial conditions are thus as follows: we have a 
static equilibrium kink and at i = the force Eq. (pn|) is 
turned " on" . 

In contrast to the majority of studies on the soliton 
ratchets we study the kink ratchets in the models that 
include no viscosity terms. We thus restrict ourselves to 
the case of driving force with a small amplitude {A < 
0.04) and with frequency lying outside the phonon band 
(more precisely, below the phonon spectrum), otherwise 
phonon modes will be excited and the analysis of kink 
ratchets will become more complicated. Study of the 
undamped ratchet makes it possible to directly measure 
the force driving the kink. 

For the models with damping terms it is customary to 
measure the efficiency of ratchet by the averaged velocity 
of steady motion of the soliton. This is not applicable 
to our case and instead, we measure the acceleration of 
the kink, a. As it will be shown, this approach works 
well for DKGM2, which is the primary subject of present 
study. This is demonstrated by the numerical results 
presented in Fig. 3] where in (a) we show two examples 
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FIG. 5: Same as in Fig. |4]but for DKGMl, where kink expe- 
riences PNp. 

of kink's trajectories (lines oscillating with the frequency 
of driving force, Vl = 0.5) and the least-square fit to these 
lines by square parabola 

x{t) = at^ + VQt + xo, (21) 

where a is the net acceleration of the kink and vq , xq are 
the kink's initial velocity and coordinate, respectively. 
One can see that kink's trajectories are fitted very well 
by the square parabola, meaning that the motion of kink 
is uniformly accelerated within the studied time domain. 

The two trajectories shown in Fig. 5(a) correspond to 
different initial phases of the driving force Eq. ([20|) . 
while all other parameters are same: A = 0.04, Q. = 0.5, 
h — 0.6. In the case of = kink does not get initial 
momentum (in this case vq — 0) while in the case of 
</5 = O.Itt it does (fo ^ 0). However, acceleration a of the 
kink in both cases is nearly same. In the panels (b) and 
(c) of Fig. 5 we plot the acceleration a and the initial 
velocity vq of the kink, respectively, as the functions of 
the initial phase of driving force, ^p. It can be seen that 
fo changes noticeably with tp but its average over is 
zero. On the other hand, the acceleration of the kink is 
practically independent of ip and in the rest of the paper 
we set ip = Q. 

We have also checked how the kink's acceleration a 
depends on the initial position of the static kink with 
respect to the lattice, xq, and found that a practically 
does not depend on xq for the chains with h = 0.3, h = 
0.6, and h — 0.9. Even for the largest studied value of 
h — 0.9 the difference between a measured for kinks with 
different xq was within the numerical error. 

In Fig. [5] we show same as in Fig. 2] but for DKGMl 
where kink experiences PNp. The results in this case are 
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FIG. 6: Acceleration of kink in DKGM2 as the function of 
frequency of driving force for different values of the force am- 
plitude, A, as specified near each curve. Initial phase of the 
force </5 = and the discreteness parameter h = 0.6. Ver- 
tical solid lines show the frequencies of the kink's internal 
modes and the vertical dashed line shows the lower edge of 
the phonon band. 



strikingly different. Kink's trajectories are now irregular 
and their least-square fit by square parabola does not 
make much sense. Nevertheless, in the panels (b) and 
(c) we present the values of a and vq obtained from such 
fit of trajectories corresponding to various ip. Both a 
and V{) vary irregularly. We thus conclude that presence 
of PNp largely affects the ratchet dynamics of kink in 
our settings. In contrast to the DKGM2, where PNp 
is absent, motion of kinks in DKGMl is not uniformly 
accelerated, at least for the range of parameters studied 
in this work, i.e., for rather small amplitude A of the 
driving force. This is true already at moderate degree of 
discreteness, h = 0.6, and the influence of PNp increases 
with increase in h. 

Now we turn back to the DKGM2 and study the in- 
fluence of the driving force parameters A and and the 
discreteness parameter h on the uniformly accelerated 
dynamics of the kink. 

Results presented in Fig. [S]were obtained for h = 0.6. 
Here we plot how the kink's acceleration a depends on 
the driving force frequency at different values of the 
amplitude of the force A, as indicated near each curve. 
Vertical solid lines show the frequencies of the kink's in- 
ternal modes and the vertical dashed line shows the lower 
edge of the phonon band. It is readily seen that the ac- 
celeration of the kink increases by one or even two or- 
ders of magnitude (note the logarithmic scale for the or- 
dinate) when the driving force frequency approaches 
the frequency of kink's internal mode luim — 1-32 (also 
note a smaller peak at u;im/2). This finding agrees well 
with earlier observations on the role of the kink's inter- 
nal modes on ratchet dynamics in the underdamped case 
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FIG. 7: Acceleration of kink in DKGM2 as the function of 
frequency of driving force for different values of the discrete- 
ness parameter ft, as specified in each panel. Amplitude of 
the ac driving force is ^4 = 0.04 and initial phase of the force 
is (/p = 0. Vertical solid lines show the frequencies of the kink's 
internal modes and the vertical dashed line shows the lower 
edge of the phonon band. 



P, nil, I23, m, [2^ . One can also see that a increases when 
approaches kink's another internal mode frequency 
'^IM = 1-73. However, this increase can also be attributed 
to the fact that the phonon band edge (wi = 1.795) is also 
approached. A special investigation is required to clarify 
which of these two factors plays major role in increasing 
a. Looking at Fig. [6l one can also notice that the in- 
crease of the driving force amplitude A by one order of 
magnitude has resulted in the increase in a by two orders 
of magnitude and thus, a ^ A^. Note that the scaling 
rule (v) ~ A^, where (v) is the averaged stationary kink 
velocity, has been reported for the soliton ratchets 1] . 

Finally, we discuss the influence of the discreteness pa- 
rameter h on the kink's acceleration a at various driving 
force frequencies n (see Fig. [7]). Here we set A = 0.04 
and consider the cases of relatively weak {h = 0.3), mod- 
erate {h = 0.6), and strong {h = 0.9) discreteness. Again, 
vertical solid lines show the frequencies of the kink's in- 
ternal modes and the vertical dashed lines show the lower 
edge of the phonon band. Remarkably, for fl < 1.2 the 
results are very close for all three values of h. This can 
be interpreted in such a way that for the discrete model 
supporting PNp-free kinks the ratchet dynamics is more 
like in the continuum case. The difference in the results 
that appears for > 1.2 is related to the kink's inter- 
nal modes, whose frequencies are /i-dependent [see Fig. 
in^b)]. We note that there is no numerical data in Fig. 
^c) for n > 1.325. Above this frequency there is a mixed 



influence of the two kink's internal modes and motion of 
the kink becomes different from uniformly accelerated so 
that one cannot assign any particular value of a to it. 



V. CONCLUSIONS 

Undamped ratchet dynamics of discrete Klein-Gordon 
kinks free of the Peirls-Nabarro potential was investi- 
gated numerically. For this purpose a lattice with asym- 
metric on-site potential was constructed. 

It was found that, typically, in the presence of single- 
harmonic AC driving and in the absence of damping, 
PNp-free kink dynamics is uniformly accelerated until 
its velocity becomes too large and radiation losses start 
to contribute to the dynamics. 

Our main finding is that discrete kink ratchets in the 
absence of PNp, at least for relatively small amplitude 
of driving force, is very much different from the conven- 
tional discrete kink ratchets experiencing PNp (compare 
the results plotted in Fig. 2] and Fig. [5]). 

Particularly, the acceleration of the PNp-frce kink due 
to AC driving practically does not depend on h in the 
non-resonance range of the driving frequency Q (see Fig. 
[7]). Indeed, in the range of 0.1 < fl < 1.2 we have practi- 
cally same a{fl) dependence for relatively weak {h = 0.3), 
moderate {h = 0.6), and strong {h = 0.9) discreteness. 
Is it is well-known that typically, physical properties of 
a discrete system are extremely sensitive to the discrete- 
ness parameter h. This unusual result can be explained 
by the fact that for the static kinks in DKGM2 PNp is 
precisely equal to zero. 

Influence of h on a{^l) appears only through the influ- 
ence of the kink's internal modes whose frequencies are 
/i-dependent [see Fig. EJb)]. 

We also confirm earher findings [H, HH HE HI, [ll that 
the efficiency of ratchets, measured in our undamped 
case by the acceleration of kink, considerably increases 
when the driving force frequency approaches kink's in- 
ternal mode frequency (see Fig. [HI and Fig. [7]) . We also 
observed that a ^ A^, where a is kink's acceleration and 
A is the driving force amplitude. This is similar to the 
scaling rule (v) ^ A^ reported earlier for the averaged 
soliton velocity, (v), [1]. 

Main reason for the striking difference in the kink 
ratchet dynamics observed for PNp-free and ordinary 
kinks lies in the fact that the static PNp-free kinks are 
not trapped by the lattice and they possess the zero- 
frequency translational Goldstone mode for any degree 
of discreteness [see Fig. El^b)]. 

Many interesting problems are left out of the scope of 
the present work, such as influence of damping, the effect 
of asymmetric, e.g., biharmonic AC driving, stochastic 
driving, etc. We plan to address these issues in forth- 
coming publications. 
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